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We study the process 77^J/i/)cj using a data sample of 519.2 fb^^ recorded by the BABAR 
detector at SLAC at the PEP-II asymmetric-energy e+e~ collider at center-of-mass energies near 
the T{nS) (n = 2, 3, 4) resonances. We confirm the existence of the charmonium-like resonance 
X(3915) decaying to J/'i/'w with a significance of 7.6 standard deviations, including systematic 
uncertainties, and measure its mass (3919.4 ± 2.2 ± 1.6) MeV/c^ and width (13±6±3) MeV, where 
the first uncertainty is statistical and the second systematic. A spin-parity analysis supports the 
assignment = 0"*" and therefore the identification of the signal as due to the Xco(2P) resonance. 
In this hypothesis we determine the product between the two-photon width and the final state 
branching fraction to be (52 ± 10 ± 3) eV . 
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I. INTRODUCTION 



In the last several years many new charmonium-like 
states have been observed in the mass region between 
3.7 and 5.0 GeV/c^, above the DD threshold, with prop- 
erties that disfavor their interpretation as conventional 
charmonium mesons The A:(3915) resonance, de- 

caying to the J/^uj final state, was first observed by 
the Belle Collaboration in two-photon collisions ^Q] . An- 
other resonance, dubbed F(3940), has been observed in 
the B^J/ipojK process [J, [a, The mass measure- 
ment for the y(3940) dlil, 01 is consistent with that 
of the X{3915) Q. Thus, the same particle, with a 
mass of about 3915 MeV/c^, may have been observed 
in two distinct production processes. The Z{3930) reso- 
nance has been discovered in the jj-^DD process 0, Q . 
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Its interpretation as the Xc2(2P), the first radial ex- 
citation of the charmonium ground state, is com- 
monly accepted Interpretation of the X(3915) as the 
Xco(2P) P or Xc2(2P) state fTo| has been suggested. The 
latter implies that the X{3915) and Z(3930) are the same 
particle, observed in different decay modes. However, the 
product of the two-photon width times the decay branch- 
ing fraction B for the X(3915) reported by Belle Q is un- 
expectedly large compared to other excited cc states Q . 
Interpretation of the X(3915) in the framework of molec- 
ular models has also been proposed [Tlj . 

Despite the many measurements available @, the na- 
ture of the X(3872) state, which was first observed by 
Belle [13] , is still unclear [l3| ■ The observation of its de- 
cay into jJ/ip [l3l ensures that this particle has positive 
C-parity. The spin analysis performed by CDF on the 
decay X{3872)^J/iPtt+tt- concludes that only = 1+ 
and — 2" are consistent with data flE'l. Similarly, 
a recent spin analysis performed by Belle [16] concludes 
that — 1^ describes the data as does = 2^ with 
one free parameter. An analysis of the Tr+Tr^Tr*^ mass 
distribution in the X{3872)^J/ipuj decay performed by 
BaBAR favors the spin-parity assignment = 2~ [3], 
but a J-^ = 1"*" spin assignment is not ruled out. If 
= 2^, the production of the X(3872) in two-photon 
collisions would be allowed. 

In this paper we search for the Ar(3915) and 
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X{3872) resonances in the two-photon process 
e+e"— 7'e"'"e~77— >e+e^ where J/tp^£~^£~, {£ — 
e or /i) and uj-^tt^tt^tt'^. Two- photon events where 
the interacting photons are not quasi-real are strongly 
suppressed in this analysis by the selection criteria 
described below. This implies that the allowed J^*^ 
values of any produced resonances are 0=*=+, 2=*=+, 4=*=+, 
S"*"*", 5"*"^, ... [l7| . Angular momentum conservation, 
parity conservation, and charge conjugation invariance 
then imply that these quantum numbers also apply to 
the final state. 

This paper is organized as follows. In Sec. II we give 
a brief description of the BaBAR detector. Section III is 
devoted to the event reconstruction and data selection. 
In Sec. IV we present the study of the J/i/joj system while 
in Sec. V we perform an angular analysis of X{3915). 
The study of systematic uncertainties is described in Sec. 
VI. In Sec. VII we summarize the results. 



II. THE BABAR DETECTOR 

The results presented here are based on data collected 
with the BaBAR detector at the PEP-II asymmetric- 
energy e+e" collider located at the SLAC National Ac- 
celerator Laboratory, and correspond to an integrated 
luminosity of 519.2 fb"""^ recorded at center-of-mass ener- 
gies near the T{nS) {n = 2, 3, 4) resonances. The BaBAR 
detector is described in detail elsewhere [l8|. Charged 
particles are detected, and their momenta are measured, 
by a five-layer double-sided microstrip detector and a 40- 
layer drift chamber, both operating in the 1.5 T mag- 
netic field of a superconducting solenoid. Photons and 
electrons are identified in a CsI(Tl) crystal electromag- 
netic calorimeter (EMC). Charged-particle identification 
(PID) is provided by the specific energy loss in the track- 
ing devices, and by an internally reflecting, ring-imaging 
Cherenkov detector. Muons and neutral mesons are 
detected in the instrumented flux return of the magnet. 
Monte Carlo (MC) simulated events with sample 
sizes more than 10 times larger than the corresponding 
data samples, are used to evaluate the signal efficiency 
and determine background features. Two-photon events 
are simulated using the GamGam MC generator [31]. 



III. EVENTS RECONSTRUCTION AND DATA 
SELECTION 

In this analysis we select events in which the and 
e~ beam particles are scattered at small angles and re- 
main undetected. In the 77— J-J/i/io; process, the J/ip is 
reconstructed in the i'^i^ final state, with £ — e or fi, 
while the oj is reconstructed in its dominant tt^tt^tt*^ de- 
cay mode. We only consider events where the number 
of well-measured charged tracks having a transverse mo- 
mentum greater than 0.1 GeV/c is exactly equal to four. 

Neutral pious are reconstructed through the tt''— ^77 



decay. We require the invariant mass of a 7r° candidate 
to be in the range (115-150) MeV/c^, and its energy in 
the laboratory system to be larger than 200 MeV. The 
energy in the laboratory frame of the most energetic pho- 
ton from 7r° decay is required to be smaller than 1.4 GeV 
in order to suppress tt'^'s not originating from an lo de- 
cay. We require the energy of the least energetic photon 
from tt" decay to be in the range (0.04-0.60) GeV, and 
I cosH^ol < 0.9, where 'H^o is the angle between the sig- 
nal 7r° flight direction in the laboratory frame and the 
direction of one of its daughters in its rest frame. These 
requirements are optimized by maximizing S/V S + B, 
where S is the number of MC signal events with a well- 
reconstructed 7r°, and B is the number of MC signal 
events where the 7r° is misreconstructed. The w is re- 
constructed by combining two oppositely charged tracks 
identifled as pions with one tt". The uj signal region is 
deflned as 740 < m(7r+7r-7r'^) < 820 MeV/c^. The J/^p 
is reconstructed by combining two tracks that are identi- 
fied as oppositely charged muons or electrons. The mea- 
sured electron energy is corrected to account for energy 
deposits in the EMC consistent with bremsstrahlung ra- 
diation. We require the vertex fit probability of the two 
leptons to be larger than 0.1%. The J/tp signal region 
is defined as 2.95 < m(e+e") < 3.14 GeV/cP for e+e~ 
and 3.05 < m{ii'^^~) < 3.14 GeV/c^ for ji^^i" events. 
An event with a J/ipu> candidate is constructed by fitting 
the J/ip and lu candidates to a common vertex. The 
mass is constrained to its nominal value [8| in this fit. 
Charged particles are required to originate from the in- 
teraction region. We require the vertex fit probability of 
the charmonium candidate to be larger than 0.1%. 

Background arises mainly from random combinations 
of particles from e'^e~ annihilation, other two-photon 
collisions, and initial-state radiation (ISR) processes. We 
discriminate against J/?/;7r+7r^7r° events produced via 
ISR by requiring .^^ = (p^+^. -p^oc)' > 2 (GeV/c^)^ 
where Pe+e- (Prec) is the four-momentum of the initial 
state {J/tpu final state). We define pt as the transverse 
momentum, in the e"'"e~ rest frame, of the J/ipoj candi- 
date with respect to the beam axis. Well-reconstructed 
two-photon events are expected to have a low trans- 
verse momentum pT and a small amount of EMC energy 
-E^cxtra, i-e-, energy not associated with the final state par- 
ticles. We require Pt < 0.2 GeV/c and i^cxtra < 0.3 GeV. 
Events originating from residual ISR ■0(25')— )■ J/V'tt+tt" 
decays may create fake structures in the J/ipoj mass spec- 
trum. We therefore remove events in the mass win- 
dow 3.675 < m(J/07r+7r-) < 3.700 GeV/c^, where 
TO(J/V'7r+7r^) = m{£+t-Tr+TT-)-m{i+£-) + m{J/iPf^'^ 
and m{J/ip Y'^'^ is the nominal J/ip mass |8|]. 

The J/ipuj signal region is defined as the intersection 
of the J/ip and lo signal regions defined above. In about 
10% of the events we find more than one candidate, and 
we select the one having the lowest px value. We obtain 
95 events in the Jjipio signal region. 



7 



IV. STUDY OF THE J/i/^a; SYSTEM 

Figure [1] shows the pr distribution for the selected can- 
didates, obtained by applying the above requirements 
with the exception of that on px- The distribution is 
fitted with the signal pt shape obtained from MC simu- 
lation plus a combinatorial background component, mod- 
eled using a second-order polynomial function with free 
parameters. The number of events from combinatorial 
background in the px < 0.2 GeV/c region is 4 ± 3. 




FIG. 1: The pr distribution of selected candidates (solid 
points). The solid histogram represents the result of a fit 
to the sum of the simulated signal (dashed) and background 
(dotted) contributions. 



Figure shows the distribution in the m{£~^£~)- 
m(7r+7r~7r*') plane of events that satisfy the selection cri- 
teria, except for the J/ip and uj mass selections. The fig- 
ure also shows the definitions of signal and background 
regions, indicated by the tiles labeled 1-9. The signal re- 
gions correspond to tile 5. Figures ^a) and (b) show 
m{£^£^) and m(7r~''7r~7r'^) for events in the w and J/ip 
signal regions, respectively. As a consistency check, we 
assign an w-Dalitz-plot weight 0] to events in the J/ipijJ 
signal region. The procedure makes use of the uj decay 
angular distribution. The helicity angle 9 is the angle 
between the tt"*" and directions in the tt^tt^ reference 
frame. The cos 9 distribution is proportional to sin^ 0, 
and the ut signal is projected by giving the z**^ event 
weight Wi — |(1 — 3 cos^ 9i). The sum of the w-Dalitz- 
plot weights is consistent with the number of events in 
the J/ijjuj signal region, thus consistent with the hypothe- 
sis that most of the observed events do indeed arise from 
true UJ — > Tr+TT^Tr*^ decays. 

To improve the mass resolution, we define the recon- 



structed J/tpiv mass as m{J/ipu!) — m{£^£ 



-vr") - 



m{£~^£~) + m{J/tjj)^^'~^ ■ The non- J/ipuj background is es- 
timated from the J/tp and w sidebands defined in Fig. [21 
The u! sidebands are defined as [0.55,0.59] GeV/c^and 
[1.00,1.04] GeV/c^. The J/^p sidebands are defined as 



o 

0) 




3.2 3.3 
m(l*l") (GeV/c^) 

FIG. 2: Event distribution (solid points) in the m(7r^7r~7r°) 
versus m{£^£~) plane for the (a) e"'"e~ and (b) j-t'^ fj,~ decay 
mode of the J/^p . We also show the J/^p ui signal region (tile 
5) and sidebands (tiles 1-4 and 6-9). 
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FIG. 3; Data (sohd points) and normalized MC (histogram) 
distributions of (a) m{£'^£~) for events in the oj signal region, 
and (b) m(7r^7r~7r°) for events in the J/ip signal region. 



[2.805,2.900] GeV/c2 and [3.170,3.265] GeV/c^ for the 
e+e" channel, [2.970,3.015] 



GeV/c^ and [3.170,3.215] 
GeV/c^ for the channel. With these definitions, 

each sideband size is half of the signal size. The m{J/ipLd) 
spectrum of this background in the J/ip lo signal region is 
obtained by 5(5) = B{2) + 5(4) + 5(6) -I- 5(8) - (5(1) + 
5(3) + 5(7) + B{9)), where B{i) is the m{J/iPuj) spec- 
trum in the z*^ region shown in fig. [51 The estimated 
background from this method is 5 ± 3 in good agreement 
with the estimate from the fit to the px distribution. The 
residual background from ■0(25)— >■ J/'07r"'"7r~ decay is es- 
timated by using the values of the integrated luminosity, 
MC efficiencies, the cross section for ip{2S) production 
in ISR events [20|, and the nominal branching fractions 
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for the relevant i^(2S) and J/ip decays Q. The expected 
number of background events from such process is smaUer 
than 0.9 at 90% confidence level (CL). 

The detection efficiency depends on m{J/ipuj) and 
where is the angle between the direction of the posi- 
tively charged lepton from J/-0 decay (^+) and the beam 
axis in the J/ipLu rest frame. Since we select events in 
which the and e~ beam particles are scattered at 
small angles, the two-photon axis is approximately the 
same as the beam axis. Therefore we use the beam axis 
to determine 

We parameterize the efficiency dependence with a two- 
dimensional {m{J/ipuj), 0*^) histogram. We label MC 
events where the reconstructed decay particles are suc- 
cessfully matched to the generated ones as truth-matched 
events. The detection efficiency in each histogram bin is 
defined as the ratio between the number of truth-matched 
MC events that satisfy the selection criteria and the num- 
ber of MC events that were generated for that bin. 

The m{J/ipuj) spectrum is shown in Fig. 21 where 
each event is weighted to account for detector efficiency, 
which is almost uniform as a function of the J/^juj mass. 
The event weight is equal to e/e{m{J/ipuj),6i), where 
e{m{J/tpLd),9g) is the m{J/ipuj)- and 6'|-dependent effi- 
ciency value and e is a common scaling factor that en- 
sures all the weights are 0(1), since weights far from 
one can cause the estimate of the statistical uncertainty 
to be incorrect [2lj. We observe a prominent peak near 
3915 MeV/c^ over a small background. No evident struc- 
ture is observed around 3872 MeV/c^. 

We perform an extended unbinned maximum- 
likelihood fit to the efficiency-corrected m{J/ipu!) spec- 
trum to extract the resonance yield and parameters. In 
the likelihood function C there are two components: one 
for the X{3915) signal and one for the non-resonant 
J/ipuj contribution (NR). The probability density func- 
tion (PDF) for the signal component is defined by the 
convolution of an 5-wave relativistic Breit-Wigner dis- 
tribution with a detector resolution function. The NR 
contribution is taken to be proportional to Vbgim) — 
p*{m) X exp[—Sp*{m)], where p*(rn) is the J/ip momen- 
tum in the rest frame of a J/tp ui system with an invariant 
mass m, ^ is a fit parameter, and m = m{J/^uj). The 
signal and NR yields, the X(3915) mass and width, and 
5 are free parameters in the fit. 

We use truth-matched MC events to determine the 
signal PDF detector resolution function. The signal 
detector-resolution PDF is described by the sum of two 
Gaussian shapes for the X(3915) and the sum of a Gaus- 
sian plus a Crystal Ball function [22] for the X(3872). 
The parameters of the resolution functions are deter- 
mined from fits to truth-matched MC events. The 
widths of the Gaussian core components are 5.7 MeV and 
4.5 MeV, respectively, for X(3915) and X(3872). No sig- 
nificant difference in the resolution function parameters 
is observed for the different J/-0 decay modes. The pa- 
rameters of the resolution functions arc fixed to their MC 
values in the maximum- likelihood fit. 



The fitted distribution from the maximum-likelihood 
fit to the efficiency-corrected m{J/'ipuj) spectrum is shown 
in Fig.m We observe 59 ±10 signal events; the measured 
X(3915) mass and width are (3919.4 ± 2.2) MeV/c^ and 
(13 ± 6) MeV, respectively, where the uncertainties are 
statistical only. We add an X{3872) component, mod- 
eled as a P-wave relativistic Breit-Wigner with mass 
3872 MeV/c^ and width 2 MeV S], convoluted with 
the detector resolution function. No significant change 
in the result is observed with the addition of this com- 
ponent, whose yield is estimated to be 1 ± 4 events. 
An excess of events over the fitted NR is observed at 
m{J/ipuj) ~ 4025 MeV/c^. If we add a resonant compo- 
nent in the likelihood function to fit this excess, modeled 
as a Gaussian having free parameters, we obtain a signal 
yield of 5 ± 3 events. 



i 15 
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FIG. 4; The efSciency-corrected m(J/7/)aj) distribution of se- 
lected events (solid points). The solid line represents the 
total fit function. The dashed line is the NR contribution. 
The shaded histogram is the non-J/ijjui background defined 
in the text as 5(5) and estimated from sidebands. The verti- 
cal dashed (red) line is placed at m(J/i/) cj) = 3.872 GeV /(? . 



V. ANGULAR ANALYSIS OF THE X(3915) 

We first attempt to discriminate between — 0* and 
= 2"*" by using the Rosner predictions. In addi- 
tion to the previously defined 0*^ we consider the follow- 
ing two angles: 0* defined as the angle between the nor- 
mal to the decay plane of the lo (n) and the two-photon 
axis, and 9in defined as the angle between the lepton 
from J/ip decay and the uj decay normal (see Fig. [5]). 
To obtain the normal to the uj decay plane we boost the 
two pious from the uj decay into the uj rest frame and 
obtain n by the cross product vector of the two charged 
pious. A projection of the efficiency values over cos0^ in 
the X(3915) signal region is shown in Fig.[SKa). The pro- 
jections of the efficiency over the angles 6** and 6'/„ are 
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FIG. 5: Diagram illustrating the reference frames involved in 
the definition of angular variables. 



shown in Figs. [H^b) and Fig. [He). The efficiency distri- 
butions are not uniform and are parameterized by fifth- 
order polynomials. The cos6'| cos0*, and cos 9in distri- 
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FIG. 7: The efficiency-corrected distributions of selected 
events in the X(39f5) signal region 3890 < m(J/i/)cj) < 
3950 MeV/c^ (solid points), (a) cosdl (b) cos 6^, (c) cos 6i„, 
and (d) cos6h- The solid (red) line represents the expected 
distribution for the = 0^ assignment and the dashed 
(blue) line for the = 2+. 




FIG. 6: The efficiency distributions in the X(39f5) signal 
region 3890 < m{J/ipLu) < 3950 MeV/c^ (solid points) as 
functions of: (a) cosO^ (b) cos 6^, (c) cos6';„, (d) cos9h, and 
(c) . The curves show the results from the fits described in 
the text. 



butions are sensitive to the spin-parity of the resonance. 
We assume that for = 2+ the dominant amplitude 
has helicity 2. This is i n ag reement with previous char- 
monium measurements [2J-|26[, and theoretical predic- 
tions [53, [15 . The expected functional forms under this 
hypothesis are summarized in Table ID Figures [71[a),(b), 
and (c) show the efficiency-corrected cos6'| cos0*, and 



cosOin distributions for events in the X(3915) signal re- 
gion, defined by 3890 < m{J/^w) < 3950 MeV/c^. Since 
the background is small, we assume that all the events 
come from X(3915) decay. The distributions for data 
are compared with the expected curves for = 0^ and 
= 2+. The resulting for each distribution is re- 
ported in Table m In all cases the = 0^ expectations 
describe the data better than the — 2+ ones and this 
is particularly true for the cos 0* distribution. In the lat- 
ter case probabilities for = 0^ and = 2+ are 
respectively 64.7% and 9.6 x 10^^% respectively. We con- 
clude that the data largely prefer = 0* over = 2+. 

The spin-0 hypothesis can be further tested by exam- 
ining the cos 9h distribution, where 9^ is the angle formed 
by the J/^ momentum in the J/tpuj rest frame with re- 
spect to the J/i}}UJ direction in the laboratory frame. The 
efficiency distribution as a function of cos 9^ is shown in 
Fig.lHd), where it is parameterized by a third-order poly- 
nomial. The cos 0ft, distribution in the X(3915) signal 
region, corrected for efficiency, is shown in fig. El^d) and 
is compared with the uniform distribution expected for 
the spin-0 hypothesis. The resulting x^/NDF is 12.2/9 
and we conclude that this test also supports the spin-0 
assignment. 

We attempt to discriminate between = 0^ and 
J'^ = 0"*". For this purpose, we define the angles, 6'„, 9i, 
and To define these angles, we first boost all the 4- 
vectors into the J/i/joj rest frame. We define 6'„ to be the 
angle between the normal to the uj decay plane n and the 
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TABLE I: Functional shapes and ^or the different spin 
hypotheses. NDF=9. 



Angle 


jP = 0" 




jP ^ Q± 


= 2+ 


"I 






1 




X 






1 1 9 


1 fi Q 








1 


sin^ f?* 


X 






f\ Q 

o.y 


DO.y 


9} 

'-'In 






sin^ 9, 


7 _ cos^ 6"? 


x' 






12.5 


18.0 


Oh 






1 




x' 






12.2 




dn 


sin*^ On 


1 






x' 


77.6 


16.3 








1 + cos^ Oi 


1 






x' 


8.7 


8.3 






<t)i 2 


~ cos(2 cos (j) 


) 2 + cos(2 cos (j) 


) 




x' 


21.7 


9.6 







Lo direction in the J /tpuj rest frame. The efficiency distri- 
bution as a function of cos 0„ (not shown) is consistent 
with being uniform. 

For J/ip decay, we first boost the f"*" to the J/-0 rest 
frame. We define 9i as the angle between the and 
the direction of the J/ip in the J/'f/'w frame. The effi- 
ciency distribution as a function of cos Oi (not shown) is 
consistent with being uniform. 

Next we define a coordinate system as follows (see 
Fig. For (jj decay, we choose the z axis along the 
Lo momentum vector, and represent the w decay in terms 
of its decay plane normal n. The cross product z x n gives 
the y-axis direction. Then we define the x-axis vector by 
y y. z. The x — z plane, by construction, contains the lo 
decay plane normal. 

We now specify the J/ij} decay coordinate system in 
terms of the unit vectors defined for u decay. We define 
z' = —z, x' = —X, and tf = y so that is along the 
normal to the plane containing the normal to the decay 
plane of the lu. Next we define the J/ip decay plane nor- 
mal n' as the cross product of the ^+ in the J/i/j rest 
frame and the z' vector. By construction, n' is in the 
x' — y' plane. Then we compute the angle (pi as the angle 
between the J/ip and uj decay plane normals. 

The efficiency distribution as a function of (pi is shown 
in Fig. and is fitted using the function e(0/) = 1 — 
c • cos2(/);, where c is a free parameter. 

It can be shown that the full angular distribution for 
= 0^ can be written as: 



3 /I /I , , = It- sin^ 6'„[l+cos6';^+sin^ 6'rcos 20;] 

rfcos On acos Oidcpi 647r 




0.5 1 

C0S9, 



FIG. 8: The efficiency-corrected distributions of selected 
events in the X(3915) signal region 3890 < m(J/i/)aj) < 
3950 MeV/c^ (solid points). The (a) cos6l„, (b) cos 61;, and 
(b) distributions are compared with — 0^ (solid line, 
red) and = 0~ (dashed line, blue) expectations. 



angular distribution is given by 
dN 



dcos Ondcos 



= [2 sin^ 01 cos^ 6'„4 

sin^ 6'n ■ (1 -l-cos^e*; - 
sin 2di cos (pi sin 20„] . 



sin^ 9i cos 20;)- 



(1) 

For = 0+, assuming no D-wave, the normalized 



(2) 

Eqs. (1) and (2), when projected onto the different an- 
gles, give the functional expectations shown in Table U 
and presented in Fig. [S] The resulting for all the 
distributions are summarized in Table ID In all cases the 
= 0+ hypothesis gives a smaller than the = 0~ 
hypothesis and this is particularly true for the cos 6n dis- 
tribution. In the latter case probabilities for = 0+ 
and = 0~ are 6.1% and 4.8x10^^^% respectively. 
We conclude that the — 0+ assignment is largely pre- 
ferred over the — 0^ assignment. 

We observe no correlation between any angles consid- 
ered in this analysis except for 0; which is strongly cor- 
related with Bin- 



VI. SYSTEMATIC UNCERTAINTIES 

Several sources contribute to systematic uncertainties 
on the resonance yields and parameters. Systematic un- 
certainties due to the functional forms chosen for the 
PDF parametrizations and fixed parameters in the fit 
are estimated to be the sum in quadrature of the changes 
observed when repeating the fit varying the fixed param- 
eters by ±1 standard deviation (a). Since the X(3915) 
spin assignment is unknown, we repeat the fit by pa- 
rameterizing the X(3915) signal as the convolution of a 
P-wave relativistic Breit-Wigner with the detector reso- 
lution function. The changes in the fit results are taken as 
the systematic uncertainty. We examine the dependence 
of the fit results on the fit range, varying the boundary of 
the fit from the nominal value of 4.2 GeV/c^ (see Fig. 2]) 
to either 4.1 or 4.3 GeV/c^. We take as the system- 
atic uncertainty the largest among the observed differ- 
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ences in the fit results. The uncertainty on the absolute 
mass scale is studied by measuring the difference between 
the observed and nominal J/ili mass in a K^tt'^it^t:^ 
ISR-enriched control sample The K+ R- 11+ -r- tt^ 

final state has the same number of charged and neutral 
particles as J/ipuj. The observed difference in mass is 
(—1.1 ± 0.8) MeV/c^. We take the sum in quadrature 
of this shift with its uncertainty as a systematic uncer- 
tainty. Previous studies show that MC events have a 
better mass resolution than data [2^ . The effect of pos- 
sible data/MC differences in the m{J/ipuj) resolution is 
estimated by increasing the width of the resolution func- 
tion core component by 20%. The uncertainty due to the 
use of efficiency weights to correct the m{J/^ui) spectrum 
is estimated with simulated experiments. In each exper- 
iment, we randomly modify the efficiency weight accord- 
ing to its statistical uncertainty. We then fit the resulting 
mass spectra and plot the resulting yields and resonance 
parameters. The resulting spreads give the systematic 
uncertainties on these quantities. We find that the fit 
bias on the yield is negligible. 

The X(3915) signal significance is 7.6f7, calculated 
from —2 In (£0 / J^max), where £0 and C are the likelihoods 
of the fits with and without the resonant component, re- 
spectively. The difference in the number of degrees of 
freedom is taken into account. Systematic uncertainties 
are incorporated into the likelihood function by convolv- 
ing it with a Gaussian with mean equal to zero and width 
equal to the systematic uncertainty on the yield. 

The product between the two-photon coupling T-^^ 
and the resonance branching fraction B to the J/tjiu) fi- 
nal state is measured using 473.8 fb~^ of data collected 
near the T(4S') energy. The efficiency-weighted yields for 
the resonances, the integrated luminosity near the T(4S') 
energy, and the branching fractions B{J/ip^£~^£~) — 
(5.94±0.06)% H and B{uj^n+n-n°) = (89.2±0.7)% [8] 
are used to obtain Tj^ x B using the GamGam generator. 
In this calculation, the X(3915) parameters are fixed to 
the values obtained from the fit. 

The uncertainties on the weighted signal yield de- 
scribed above are taken into account in the Ty^y x B 
systematic error. Systematic uncertainties on the effi- 
ciency due to tracking (0.3% per track), tt" reconstruc- 
tion (3.0%) and PID (0.1% per pion, 0.8% per lepton) 
are obtained from auxiliary studies. The uncertainty on 
the luminosity is 1.1%. The uncertainty on the nominal 
J/ip and w branching fractions used in the calculation are 
propagated in the F^^ x B error. The GamGam calcula- 
tion has an uncertainty of 3% Q . 

Since no significant X{3872) signal is observed, we 
determine a Bayesian upper limit (UL) at 90% CL on 
Tryj X B, assuming a uniform prior probability distribu- 
tion. The upper limit for Tjj x B is thus computed ac- 
cording to: 

/ L{T-^^ X B)d{T.^-y xB) = 0.90 
Jo 

where L{T^~f x B) is the likelihood function for F.^-^ x B. 



For a J = resonance, the resulting value of 
F^.^[X(3915)] X B{X{3915)^J/'tpuj) is (52 ± 10 ± 3) eV 
where the first uncertainty is statistical and the second 
systematic. For completeness we also report the value 
for J = 2: (10.5 ±1.9 ±0.6) eV. For X (3872), we obtain 
r^.^[X(3872)] X B{X{3872)^J/ijoj) < 1.7 eV at 90% CL, 
assuming J — 2. 

VII. SUMMARY 

In summary, we confirm the observation of the 
charmonium-like resonance X(3915) in the jj^J/ipuj 
process, with a significance of 7.6a, including systematic 
uncertainties. The measured mass and width are 

m[X(3915)] = (3919.4 ± 2.2 ± 1.6) MeV/c^ 
r[X(3915)] = (13±6±3)MeV, 

where the first uncertainty is statistical and the sec- 
ond systematic. These measurements are consistent with 
those previously reported by Belle for the same pro- 
cess ji] and by BaBAR [| and Belle [| for B^J/i}^ujK. 
A detailed angular analysis has been performed. We find 
that the data largely prefer ~ 0^ over = 2+. 
In this hypothesis, — 0+ is largely preferred over 
= 0^ and this would identify the signal as being 
due to the Xcoi'^P) resonance. The mass of Ar(3915) is 
consistent with the result of the potential model, which 
predicts the mass of the first radial excitation Xco to 
be around 3916 MeV according to Godfrey-Isgur rela- 
tivized potential model [s^l- The product F-y-y[Ar(3915)] x 
B[X (3915)^ J/tpuj] is also measured. The value for J = 
(relatively large compared to charmonium modelpredic- 
tions) is consistent with that reported by Belle Q. This 
product, also computed in this analysis for J = 2, is 
smaller than the corresponding value obtained by Belle. 
We have also searched for the 77— >-Ar(3872)— 7>J/-0a; pro- 
cess, but no significant signal is found. 
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